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Starting from the light-cone plasma (LCP) distribution of gluons produced in proton-proton col- 
lisions, we derive the invariant charged-particle yield dN / dyd?p±_ in heavy-ion collisions. Multiple 
scattering of partons in the other nucleus leads to p±_ broadening, which depends on the central- 
ity of the collision. The resulting change of the LCP parameters determines the distributions of 
charged particles. The charged-particle multiplicity dN/drj at r] = normalized to the number 
of participating nucleons has a universal centrality dependence, common to RHIC and LHC. This 
universal increase is of geometrical origin. Such a universality can be broken by the energy depen- 
dent coherence effect, shadowing and mutual boosting of the saturation scales. However, we found 
these effects to essentially compensate each other, leaving the centrality dependence of multiplicity 
approximately universal. 



I. INTRODUCTION 

In a recent publication T a universal multiplicity dis- 
tribution for produced partons has been deduced from 
maximum-entropy theory. Given the experimentally 
measured total transverse energy E± and the sum rule 
requiring the total light-cone fractions x of partons in 
the proton to be unity, one obtains a Bose-distribution 
n{x,p±) for gluons with two Lagrange parameters gov- 
erning its X and p± dependence. The underlying picture 
is simple. In soft high energy collisions partons interact 
by exchanging small transverse momenta and become lib- 
erated subsequently. The resulting hadronization prod- 
ucts reflect the parton distributions via hadron-parton 
duality. This picture has been worked out for proton- 
proton (pp) and nucleus-nucleus (AA) collisions in ref. 
[T]. It provides a satisfactory explanation for soft par- 
ticle production in pp collisions throughout the acces- 
sible rapidity range. Nuclear multiplicity distributions 
dN/dycfp± can be calculated by scaling the pp distri- 
bution with the number of participating nucleons and 
taking into account the higher mean transverse momen- 
tum observed in experiment. In this paper we want 
to theoretically calculate this increased mean transverse 
momentum from parton multiple scattering. 

To describe the inclusive cross section of charged par- 
ticles produced in pp collisions. 
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straints. We consider the light-cone momenta of the par- 
tons with energies e and longitudinal momenta pz relative 
to the light-cone momentum of the incoming proton with 

{E,Pz,oy. 
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we use a statistical distribution function which maximizes 
the entropy of the produced partons given certain con- 



The maximum-entropy distribution, i.e., the LCP dis- 
tribution, has to satisfy the following two requirements: 
The first requirement concerns light-cone momentum 
conservation; i.e. the sum of all partonic light-cone frac- 
tions is unity. The second constraint is given by the mea- 
sured total transverse energy of the produced partons. 
For symmetric collisions in the center-of-mass (cm) sys- 
tem, we can consider the multiplicity distribution in each 
hemisphere separately. The resulting individual distribu- 
tions function have the form of a Bose-Einstein distribu- 
tion depending on transverse momentum and light-cone 
fractional momentum x. It is important to emphasize 
the light-cone property of the maximum-entropy distri- 
bution. The dynamics of collisions at high energies is 
governed by a the light-cone Hamiltonian which is boost 
invariant and determines the distribution functions de- 
pending on transverse momentum and light-cone frac- 
tions which label the eigenstates. The density matrix 
resulting after the collision can be built up from an in- 
coherent mixture of such multi-parton states. Without a 
boost invariant formulation one cannot define a number 
density of partons since in each reference system it will 
be different. In the rest system of the proton the gluons 
sit in the clouds around the constituent quarks and the 
strings holding them together, whereas in a fast-moving 
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proton the virtual gluons materialize as partons carrying 
a sizable momentum fraction. It is this physics which 
also determines the maximum-entropy state. 

In collisions the projectile and target have to interact, 
and this interaction takes place between the soft partons 
with low Bjorken x. They acquire transverse momenta; 
therefore, the produced total transverse energy repre- 
sents an important input for the maximum-entropy dis- 
tribution. Boost invariance applies to the change of the 
measured inclusive distribution in the cm system under 
a change of cm energy. Once the distribution is formu- 
lated with the variables p_L and x it can be applied to any 
cm energy. In reality, the mean transverse momentum 
will increase with increasing cm energy and thereby also 
modify the softness of the particles in x. The maximum- 
entropy principle can describe a real nonequilibrium state 
formed in the collision or an average over many collisions. 
It is probable that each individual collision has a different 
transverse momentum spectrum, but on average we ob- 
serve an effective transverse temperature characterizing 
the cross section. 

The question arises: Is such a formulation really neces- 
sary, since we have the very successful Bjorken picture [2] 
which also includes boost invariance. Indeed in the cen- 
tral region where the rapidity distribution is flat, there 
is a considerable overlap between the two formulations, 
so the effective transverse temperature represented by 
the parameter A in the above distribution is related to 
the multiplicity distribution for y = in a similar way. 
There are also differences. In the light-cone formulation 
the rapidity distribution is not flat. Furthermore, the 
applicability of hydrodynamics is a fundamental ingre- 
dient of the Bjorken picture. From a very early time 
To ^ 0.2 fm on, a locally equilibrated fluid has to be as- 
sumed which expands longitudinally. More complicated 
modern hydrodynamical calculations [3j allow also trans- 
verse expansion with viscosity and calculate the radial 
flow parameters which seem to support this picture. A 
detailed time evolution based on the thermal equation- 
of-state results. However, it has to be carefully checked 
whether the many assumptions underlying these calcula- 
tions are really reflected in the data. Minimal assump- 
tions based on conservation laws, complemented by few 
key observables may also reproduce the main features of 
the experiments. 

We think that the maximum-entropy distribution can 
serve such a purpose. Building on the dynamics of the 
underlying pp events, nuclear collisions represent a su- 
perposition of pp events. Multiparton collisions in the 
evolving larger system increase the partonic transverse 
momenta which lead to higher hadronic momenta in nu- 
clear collisions. We want to understand the microscopic 
origin of this transverse broadening and its energy de- 
pendence as a function of the centrality of the collision. 
Thereby we can separate the geometrical and dynamical 
aspects of these processes. 

The outline of the paper is as follows. In section II 
we review the features of the longitudinal LCP of glu- 



ons which we need. In section III we calculate the mean 
density profile a "beam parton" experiences in a nuclear 
collision from averaging over "target nucleons" and eval- 
uate the resulting p±^ broadening. In section IV we con- 
sider the effects of mutual boosting of saturation scales 
and shadowing in nuclear collisions and compare with 
the simple parametrization of section III. Finally, in sec- 
tion V the effect of transverse momentum broadening or 
saturation on the central pseudo-rapidity distributions 
is calculated using the relevant LCP distributions. The 
results are discussed in Section VI. 



II. THE LCP DISTRIBUTION FROM THE 
MAXIMUM-ENTROPY PRINCIPLE 

The maximum-entropy principle for the multiplicity 
distribution in pp collisions can be best derived from 
the entropy of a Bose system which is not in thermal 
equilibrium but subject to constraints. The entropy of 
the system is proportional to the logarithm of the inte- 
grated phase space. In the following we set the Boltz- 
mann constant k — \ and h — \. On the light cone, 
phase space includes the transverse spatial coordinate 
h±^ , the transverse momentum , the longitudinal light- 
cone momentum and the longitudinal spatial variable 
Z- = l/2(< — z), where t and z are time and longitudinal 
coordinate. These variables are handled like in conven- 
tional thermodynamics, i.e. the phase space multiplied by 
the gluon degeneracy factor g = 2{N'^ — 1) and divided 
by the Planck constant = {2nhY gives the number of 
available quantum states G: 
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For high energies, Feynman scaling is a good phe- 
nomenological concept; therefore, we have multiplied and 
divided this expression by P_|_ — E + to obtain the 
light-cone momentum x ~ p-\-/P-\- and the longitudi- 
nal light-cone distance p = z_P_(. which are canonically 
conjugate variables. We further make the simplifying as- 
sumption that the distribution function and consequently 
the entropy are homogeneously distributed in transverse 
space. The integration over the 6_l coordinate can then 
be executed and gives the area L\. An estimate of the 
integral of the scaled light-cone distance / ^ is more 
subtle, since it is not independent on the rest of the vari- 
ables. It has to be done separately for valence and sea 
partons. Interpolating the x dependence of the longitu- 
dinal extension of valence and sea partons we obtain the 
integral over the scaled distance; cf. ref |T]: 
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The so motivated ansatz for the phase space on the light 
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cone is crucial for all further derivations: 



2 <Pp± dx 



(6) 



Gluons are bosons; therefore, they can occupy the 
phase-space cells in multiples. The binomial of the com- 
bined number of particles and states over the number of 
states gives the number of possibilities AF^, to dis- 
tribute Nx^p^ gluons, i.e. bosons, on Gx,p^ quantum 
states. The entropy of the system is defined by the loga- 
rithm of the phase space. For large particle numbers and 
quantum states one can express the entropy in terms of 
the phase space elements and mean occupation numbers 
'^x,p±_ = Nx^p^/Gx^p^ of each quantum state. By choos- 
ing the cell sizes small we convert the sums into integrals 
over the continuum variables and then vary the entropy 
under the two constraints of light-cone momentum con- 
servation and total transverse energy: 



(27r)2 



dx 
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p± n{x,pj_) = {Ej_). 
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The variation gives the LCP distribution function 
n{x,pj_) as maximum-entropy distribution: 



n{x,p±) 



1 
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(9) 



For details we refer to ref. [T]. 

The LCP distribution together with the measure gen- 
erate a gluon rapidity distribution of the following form: 



dN 



dyd'^p± (27r) 
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This is the connection of the maximum-entropy distri- 
bution with the semi-inclusive cross section. The phe- 
nomenological description of the multiplicity distribution 
has three parameters L^^X and w. The parameter A 
plays the role of an effective transverse "temperature" . 
The "softness" w determines how small the mean x be- 
comes. With increasing center-of-mass energies, we ex- 
pect that the effective transverse temperature A and w 
increase: The collision becomes "hotter" and the particle 
distributions "softer" . The effective transverse tempera- 
ture A is calculated from the mean transverse momentum 
which is equal to the ratio of the transverse energy and 
multiplicity in one hemisphere. 

To compare with experiment we assume parton-hadron 
duality and saturate the hadronic reaction products by 
pions such that the charged hadrons make up 2/3 of the 
total multiplicity. In addition, we convert the multiplicity 
as a function of rapidity into a function of pseudorapid- 
ity using the Jacobian depending on the pion mass. In 
reference jl] we have determined the best values of L±,X 
and w for pp collisions at cm energies of \/s = 0.2 TeV 
and = 2.76 TeV; cf. Table 1. 
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0.20 


1.34 


0.183 


3.44 


0.12 


2.20 


2.76 


1.28 


0.252 


6.81 


0.23 


4.56 



TABLE I. The first four columns give the cm energy -^/s, 
the transverse size L±, the effective transverse temperature 
App and the softness Wpp of the pp light-cone distributions. 
The next columns contain the mean transverse momentum 
squared of gluons (px) (GeV^) and the rapidity distribution 
dNch/drj at = of charged particles. 



It is noteworthy that the statistical weights from the 
spins and colors of gluons correctly reproduce the mag- 
nitude of the inclusive cross sections. 

In general, one would have to split the parton content 
of the proton into quarks and gluons such that the gluon 
momentum fractions do not integrate up to unity. The 
quarks and antiquarks have to be described by Fermi 
statistics. The individual momentum fractions will then 
depend on the effective saturation scale of the collision. 
It is clear that not only the momentum distribution of the 
LCP, but also the relative amount of quarks and gluons in 
this description will be different from an equilibrated 3- 
dimensional thermal quark-gluon plasma. With increas- 
ing energy the reaction products become more and more 
gluonic. This would also show up in the particle con- 
tent after hadronization. Detailed consequences of such 
a more accurate description have to be explored theoret- 
ically and experimentally. Here in this paper we stick to 
the approximation where gluons are the only constituents 
of the proton and make up the full momentum sum rule. 



III. TRANSVERSE-MOMENTUM 
BROADENING IN NUCLEUS-NUCLEUS 
COLLISIONS 

To extend the calculation to nucleus-nucleus collisions 
we assume that the multiplicity in the A-A collision 
is proportional to the number of participating nucleons 



part ■ 



drjd^pi^ 



dNi{pj_)) 



r7i^cosh^(y) dyd^pA 



(11) 



However, this is not sufficient. It is important to take 
into account the increase of the mean transverse mo- 
mentum {p±) with the number of participants. The ini- 
tial parton distributions in the projectile nucleus will be 
broadened by the interaction with the nucleons in the 
target nucleus and vice verse. This broadening will be 
calculated now. 

Transverse-momentum broadening results from coher- 
ent rescatterings with vanishing small momentum trans- 
fers. Averaging the differential cross section with the ac- 
quired transverse momentum yields the transport coeffi- 
cient which is defined as the mean transverse momentum 
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squared times the cross section, i.e. {(jp^)- This form 
has been tested in electron-nucleus and proton-nucleus 
collisions [1]. The resulting mean Ap^ is also known as 
saturation scale [51 . It can be derived from the dipole 
nucleon cross section as follows. 

In the eikonal approximation, the ejected high- 
momentum parton moves on a classical trajectory with 
impact parameter b and picks up a non-abelian phase 
factor V{b) in the background gauge field generated by 
the nucleon: 



V{b) = V exp 



15 



+ 00 



(12) 



Here V{b) is the Wilson line of the parton with impact 
parameter b relative to the proton. We use the notation 
Af^ = ^J^t", where f's are the generators of the group 
S\J{Nc) in the fundamental representation. The differ- 
ential cross section to produce a parton with transverse 
momentum p± is given by projecting the eikonal phase 
onto p± and by taking the modulus of the amplitude in- 
tegrated over all possible impact parameters 



da 
d^px 



(13) 
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Tr 



vHb')v(b) 



Hence, a fake dipole of size f± — b — b' is constructed 
from the ejected parton in the ^-amplitude and in the 
1/^-amplitude. Their trajectories are displaced from each 
other by the distance r±. The expectation values of the 
Wilson lines have to be evaluated with respect to the 
target ground state. In the dipole model, the total cross 
section for the interaction of a dipole of size fj_ with a 
target nucleon is given by 



o-dN{r±) 



d^bll- 



Tr 



0(6 + rl) 1/(6)1 \) . (14) 



We define the quantity (crpY) the integral over trans- 
verse momentum d^p± of the differential cross section 
multiplied by p^. One sees that {<jp\) is related to the 
dipole nucleon cross section: 
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This expression confirms the result derived in [6j. To 
obtain the transverse-momentum broadening of gluons 
one has to consider the dipole-nucleon cross section as a 



function of dipole size r± and dipole-nucleon cm energy 

Vs. 

Let us first consider (33)-dipoles. In the parametriza- 
tion of ref. [7] the low energy dependence of the cross 
section is parametrized via cro(s). The fast rise of the 
dipole cross section for small dipoles at high energies is 
included in the energy dependence of ^0(5): 



o-dJv(n_) = o'o(s) 



1 — exp — 



with 



ao{s) ^ 23.6 



3 r^oiS) 
8 0.44 fm" 
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mb (17) 



and 



ro(J) =0.88(s/so)~°-'' fm. 
So = 1000 GcV^. 



(18) 
(19) 



For gluons the transport-coefficient is modified by the 
color factor 9/4: 
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In pp collisions gluons are liberated with their in- 
trinsic momenta described by the light-cone distribution 
n{x,p±). In nuclear collisions this distribution is mod- 
ified. We will calculate the effect of multiple scattering 
on the transverse temperature Xaa and then insert this 
new transverse temperature in the universal light-cone 
distribution to determine the central rapidity distribu- 
tion. The effective cm energy s for multiple scattering 
has to be calculated from the collision of a transverse 
gluon with {E± ,p±,0) colliding with a proton in the op- 
posite nucleus with four-momentum (-v/s/2, 0, \/s/2): 



P±^ 



(22) 



In the following we denote the averaging over s with 
n{x,p±) by a bar. For n{x,p±) we use the parameters 
Xpp and Wpp of the corresponding light-cone gluon dis- 
tributions from Table 1 in section II. These are the pa- 
rameters which describe the multiplicity distributions for 
RHIC and LHC energies as shown in ref. PJ . We obtain 



'E{<^Plis))gn{x,pA_) 
J2n{x,p±) 

i 11.26 for = 0.2TeV 
1 21.15 for 2.76TeV. 



(23) 
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In the second step we multiply the mean transverse 
momentum in a collision with the average profile func- 
tion. This gives the mean Ap^ of gluons: 



Api = {apl)gTB{b). 



(25) 
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1. The average profile functions r(6) (see Eq. 291 for 
Au+Au at ysjvjv = 0.2 TeV (daslied curve) and Pb+Pb col- 
lisions at 'Jsnn = 2.76 TeV are shown as functions of the 
number of participants. The number of participants increases 
for more central collisions. The functional dependence of the 
resulting mean profile functions varies because the relation be- 
tween impact parameter and number of participants depends 
on cm energy. 



We calculate the averaged profile function tb(6) now: 
The profile functions 1^(6) and Tb(6) of the colliding 
nuclei A and B are obtained from the nuclear densities 
with Woods-Saxon distributions: 

TA{h) = J pA{^b^+z^)dz, (26) 
Po 



PA{r) 



1 + exp[(r - R)/a] ' 



(27) 



We use R = 6.62 (6.38) fm and a — 0.54 fm to describe 
the Pb (Au) nucleus, respectively [8 . 

The number of participants (for a definition see e.g. 
ref. [9]) depends on the inelastic nucleon-nucleon cross 
section, for which we take the values cr™^ (0.2 TeV) — 
42 mb and aj^^ (2. 76 TeV) = 64 mb from ref. pTO]. In 
AB collisions, the total number of participants, A'part(^), 
equals the number of participants in nucleus A and in 
nucleus B which explains the factor 2 in the equation 
below for symmetric AA collisions: 



^part(&) = 2 y d2^TA(s)(l-exp[-TB(6-s>K?^]). (28) 

The averaged profile function tb (b) of nucleus B is then 
obtained by using the number of participant nucleons in 
A as weight function: 



TBib) 



J d^sTAim - exp[-TB(6 - s)a'§j^])TBib- s) 

Jd^sTAim-eM-TB{b-s)a'^N]) 

(29) 

Both the number of participants and the average pro- 
file function are functions of the nucleus-nucleus impact 
parameter b. In experiment, the number of participants 
is used as a measure of the centrality of the collision. 



Therefore we eliminate the impact parameter b and rep- 
resent the mean profile as a function of the number of 
participants, i.e. centrality. The dependence of the mean 
profile function as a function of the number of partici- 
pants is rather well described by a power law: 



part ) 



part • 



(30) 



The corresponding values for the parameters are hq ~ 
0.1155 fm~^ a = 0.435 for A = 197, and no = 
0.0775 fm"^ a = 0.491 for A = 208. A plot of the numer- 
ical profiles for both nuclear collisions is shown in Fig. [T] 

For an estimate of the mean path length traversed by 
a parton in the AA collision one has to divide this profile 
density by the mean nuclear density po ~ 0.17fm~'^. One 
obtains a mean path length of approximately 10 fm for 
central collisions. Owing to the larger inelastic pp cross 
section at LHC the b distribution of the number of par- 
ticipants extends to larger impact parameters for LHC. 
Therefore, the weight of shorter trajectories is increased 
at LHC. This explains the difference of the mean profile 
functions in Fig. [ij 

Having described the relevant steps we can now cal- 
culate the transverse momentum broadening for gluons 
in nucleus- nucleus collisions for the two different energies 
y/s = 0.2 TeV and = 2.76 TeV. Neglecting the small 
differences of the average profile functions between Au- 
Au and Pb-Pb at the different energies we see that the 
transport coefficient determines the energy dependence 
of momentum broadening. This coefficient increases by 
a factor of two between RHIC and LHC, therefore the 
mean transverse momentum broadening increases by the 
same factor. In Fig. [2] and Fig. [3] the thick (blue) lines 
represent the mean value of Apj^ of gluons as a func- 
tion of the number of participants in the AA collision for 
RHIC and LHC energies calculated using the dipole scat- 
tering cross section. The two curves have almost identical 
shapes because the profile functions in Au-Au and Pb- 
Pb collisions are nearly equal; cf. Fig. 1. The shape of 
the curves is determined by the geometry of the nuclear 
collisions. 



IV. MUTUAL BOOSTING OF SATURATION 
SCALES AND SHADOWING 

In proton nucleus scattering the transverse momentum 
of a gluon in the proton propagating through the target 
nucleus B increases by the value Ap]_ , which depends on 
the profile tb(6) of the target. The resulting momentum 
broadening also defines the saturation scale Q^g of the 
intrinsic gluon distribution of the target nucleus B jlll : 



QIb - ^Pa 



(31) 



This picture of broadening strictly applies only when 
we consider nucleus-nucleus collisions as a superposition 
of proton- nucleus collisions. In nuclear collisions, how- 
ever, multiple interactions enhance the higher Fock com- 
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FIG. 2. The average broadening Ap5_ of a gluon for Au+Au 
at — 0.2 TeV is shown as functions of the number of partic- 
ipants. The number of participants increases for more central 
colhsions. The solid blue line represents the result from the 
multiple scattering calculation. The dashed red line the result 
from mutual boosting and shadowing. 
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FIG. 3. The average broadening Ap5_ of a gluon for Pb+Pb 
at ^/s — 2.76 TeV is shown as functions of the number of 
participants. The number of participants increases for more 
central collisions. The solid line represents the result from 
the multiple scattering calculation. The dashed red line the 
result from mutual boosting and shadowing. 



ponents in the wave functions of the participating nu- 
cleons, thereby increasing the saturation scale compared 
with NA collisions [T^]. The saturation scales both in 
the target and projectile nucleus are boosted to higher 
values. Intuitively, a nucleon Fock state distribution in 
the nucleus B is affected by multiple interactions, which 
enhances the contribution of multi-parton components. 
Thus, the gluon density in the bound nucleons in B in- 
creases at small x, which in turn leads to a more intensive 
interactions of these nucleons with partons from the nu- 
cleus A, i.e. enhances broadening of partons propagating 
through B. This effect of mutual boosting of saturation 
scales has been described first in [1^ and subsequently 
discussed in ref . [121 El] ■ On the other hand, multiple 
interactions at high energies are subject to shadowing, 



which works in the opposite direction, suppressing the 
magnitude of broadening and saturation scales pjj . 

For high energies the dipole-nucleon interaction is 
dominantly given by two gluon exchange. The transport 
parameter derived in the previous section can then be re- 
lated to the gluon structure function at a certain Bjorken 
X and scale Qq. In order to make this perturbative pic- 
ture coincide with the nonperturbative model based on 
the dipole cross section of section III , we determine the 
infrared scale in such a way that the two-gluon exchange 
picture on the proton reproduces the nonperturbative re- 
sult fni, 



(32) 



with 



Ql 



(33) 



Notice that both sides of the above equation depend 
on the transverse momentum of the gluon after its liber- 
ation on mass shell. We average the left and right hand 
sides of this equation over all possible pj_ values using the 
light-cone distribution npp{x,p±) generated in pp colli- 
sions with the parameters of Table. 1 to be consistent with 
Eq. (24). For scale setting we need a NLO-distribution. 



Most distribution functions like the MSTW- distribution 
|15j in NLO-order are defined for > 1 GeV^ and do 
not allow low enough virtualities to set the infrared scale 
Qq. The NLO-GRV parametrization [iB], however, ex- 
tends to low virtualities « 0.4 GeV^ and low x. We 
can use this distribution to determine the infrared scale 
Qq. We find the following solutions of the selfconsistency 
equations for RHIC and LHC energies: 



(0.2 TeV) 
g(2.76TeV) = 



= 0.35 GeV^ 
0.41 GeV^ 



(34) 
(35) 



The effective temperature and the softness parame- 
ters given in Table 1 influence the infrared scale. It is 
therefore not unreasonable that the infrared scales ob- 
tained in this way are close to the p^-scales of the mul- 
tiplicity distributions at the same energies; cf. Table 1. 
Our values for the infrared scales may in principle dif- 
fer from the pp saturation scale associated with the color 
glass parameterizations; cf. ref. |T7j. We obtain a slightly 
smaller infrared scale Qq = 0.35 GeV^ than the satura- 
tion scale Ql = 0.462 GeV^ at RHIC energy Obviously, 
the sizes of the infrared scales only marginally justify the 
use of perturbative QCD. It may be useful to consider the 
saturation of the nonintegrated gluon structure function 
in a nonperturbative formulation (cf. ref. [ISJ) based on 
string-string interactions. We will reserve such a work to 
further studies. 

The boosting effect alone overestimates the magnitude 
of broadening. It must be considered together with the 
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effect of nuclear shadowing which reduces the nuclear 
gluon distribution compared with the pure sum of nu- 
cleonic gluon distributions at small x. The nuclear gluon 
distribution does not linearly increase with the nuclear 
profile function. A simple parametrization of the results 
in Fig. 2 of reference [11] has the following form: 



gN{x)TB{h){l 

9N{x)TB{h). 



.An[TB{b)IM) 



(36) 
(37) 



The parameter a is only weakly dependent on energy. 
In the above reference the authors find for RHIC a ~ 0.15 
and for LHC a w 0.17 with tq = 0.1 fm^'^. This reduc- 
tion enters the differential equations for the the satura- 
tion scale and broadening. The increase dQ'^g of the 
saturation scale in nucleus B given by the pj_ broadening 
of a parton of nucleus A going through nucleus B rises 
with the increase of the modified profile function cItb (b) . 
A similar equation holds for dQl^: 



dQ: 



_ 37r2 
dfB{h) ~ 2 
dQlA _ 37r2 



dTA{b) 



^ - as{QlA + Ql)XA9N{xA.QlA + Ql) (38) 



+ Ql)xBgN{xB,QlB + Qo)-(39) 



Since the saturation scale of nucleus B increases 
with the saturation scale of nucleus A, the authors 
of ref. [12] have called this phenomenon mutual boost- 
ing. With increasing saturation scale also the associ- 
ated Bjorken variable xa in the gluon structure func- 
tions changes {xb is defined in the same way with index 
A^B): 



XA = 



Ql + Q 



sA 



(40) 



The right hand sides of the differential equations de- 
pend on the X values of the gluon structure functions. 
Since the structure functions become larger at smaller 
X, the increase of broadening is linked with increasing 
p±^ of the incoming gluon in the same way as the dipole 
cross where increasing means higher dipole energy and 
larger transport parameter. 

These differential equations have to be calculated with 



the initial conditions that 



QIb = for TA{b) 



TB{b) = tq (cf. Eq. 36), which guarantee that the results 
for peripheral nuclear collisions coincide with the results 
of section III obtained without boosting. We prefer these 
differential equations to the equivalent bootstrap equa- 
tions of ref. [H], because they allow to include the av- 
eraging in an easier way. In the case of symmetric AA 
collisions both equations coincide, because rs(6) = t^(6) 
and QIa = QIb- ^'^S- [s] and Fig. [sjthe dashed lines 
show the dependence of the perturbative saturation scale 
Q^A centrality expressed by the number of partici- 
pants. Comparing these values with the corresponding 
values for Ap^ from the calculation in section III one 
sees that for RHIC and LHC the mutual boosting of 
saturation scales gives approximately the same results. 



Shadowing is very important to control the amount of 
mutual boosting, which otherwise would overshoot the 
calculation from multiple scattering by a factor of two. 



V. dN/d-n DISTRIBUTIONS IN AA COLLISIONS 

Having a theoretical calculation of gluon p± broaden- 
ing or equivalently of the saturation scales in nucleus- 
nucleus collisions we can calculate the change of the 
transverse temperature A and the softness as a function 
of centrality. Because of the well defined shape of the 
LCP distributions we can relate the calculated mean p± 
values in AA collisions to dN/drj distributions. 

To describe the momentum broadening of pions one 
has to convert the information about gluons obtained in 
sections HI and IV into pion spectra. Since the transverse 
momentum broadening of gluons described in the nuclear 
multiple scattering picture is intrinsically higher than for 
pions, one cannot use here parton-hadron duality. To 
take into account the fragmentation of the gluons into 
pions, we use the fragmentation function of ref. [19] at the 
starting scale of Qq = 2GeV^. The mean relates 
the pion momentum broadening to the gluon broadening: 



D^/g{z)^mz\l~zf 



.82 



/7r/g 



= 0.11 



We use the result of Eq. ( [25] ) of section HI: 



2 

_L,7r 



{(^p'i)a'rB{b){z )^/g. 



(41) 
(42) 

(43) 



Mutual boosting explained in section IV gives saturation 
scales which have been identified with gluon momentum 
broadening. The saturation scales QIa obtained from 
the solution of the differential equation Eqs. (38l-(39) 



can also be converted into pion momentum broadening: 



(44) 



The total mean pion transverse momentum is given by 
the sum of the intrinsic pp mean transverse momentum 
{p\) given in Table 1 and the acquired (Ap^)^ above: 



(pi,tot) = (pi) + Apl^. 



(45) 



The universal form of the LCP distributions has been 
shown in ref. [T] to describe well the rapidity and p±^ 
distributions in AA collisions. Since the shape of the dis- 
tribution is fixed by the maximum-entropy principle one 
only has to modify the effective transverse temperature 
and softness parameters of the pp light-cone distributions 
to obtain the nuclear light-cone distributions. Averaged 
over all rapidities the mean transverse momenta of glu- 
ons are related to the effective transverse temperatures 
as in Table 1: 



'_L,tot) 



3.6A 



AA- 



(46) 



Using the total mean transverse momenta one obtains 
the respective values for Xaa and waa for each value of 



0. 



50 100 150 200 250 300 350 



part 



FIG. 4. The multiplicity dNch/dri at 77 = for Au+Au at 
^/s = 0.2 TeV is shown as a function of the number of partic- 
ipants. The number of participants increases for more central 
collisions. The solid line represents the result from the multi- 
ple scattering calculation. The dashed line is the result from 
mutual boosting and shadowing. Data points represent an 
average of results from the BRAHMS, PHENIX, PHOBOS, 
and STAR experiments |10| . 




50 100 150 200 250 300 350 



part 



FIG. 5. The multiphcity dN^^/dri at 77 = for Pb+Pb at 
^/s — 2.76 TeV is shown as a function of the number of par- 
ticipants. The number of participants increases for more cen- 
tral collisions. The solid line represents the result from the 
multiple scattering calculation. The dashed line is the result 
from mutual boosting and shadowing. Data points were taken 
from [lO]. 



the profile function which in turn specify the inclusive 
cross sections: 



ml dN{XAA,WAA) 



d-qd?p±_ 2 3 y m\ cosh^ y dyd?p±_ 



with 



(47) 



dN _ gL\^ 1 
dyd^px_ ^ (27r)2 exp (to_l(1/Aaa + w AAe\y\ / ^/^)) - 1 ' 

(48) 




50 100 150 200 250 300 350 



part 



FIG. 6. Relative increase factor R = {dN^+-^ /dri/{Np^rt/2)) 
I {dN^^^ / drf) of the charged-particle multiplicity in A A colli- 
sions with respect to the multiplicity in pp collisions as a func- 
tion of A'^part- The universal increase of the data (gray^RHIC, 
black=LHC) is reasonably well reproduced by the correspond- 
ing calculations (solid lines: multiple scattering only, dashed 
lines: mutual boosting and shadowing). 



In Fig. [4] and Fig. [5] we show the resulting normalized 
pseudo-rapidity distributions at ?/ = as functions of 
the number of participants. Recall that the number of 
participants has been related to the average profile func- 
tion. Figure [6] shows that the universal, i.e., ^/snn in- 



dependent, relative increase of the measured charged- 
particle muhiplicity {dN^^ /d-q/{Np^^t/'^)) with cen- 
trality is reasonably well reproduced by the calculautions. 



VI. DISCUSSION OF THE RESULTS 

We found in ref. [T] that pp and AA collisions can be 
described by light-cone distributions which are derived 
from the maximum-entropy principle and contain as pa- 
rameters a mean transverse temperature and a softness. 
Our picture underlying this description is very simple: 
After an exchange of small transverse momenta the par- 
ton content of the proton is liberated, and parton-hadron 
duality relates the produced pions to this light-cone dis- 
tribution. One can generalize this picture to AA colli- 
sions letting the semi-inclusive cross sections scale with 
the number of participant nucleons; cf. Eq. |47[ How- 
ever, this was not sufficient to explain the experimental 
data. In addition, it was necessary to use the increased 
mean transverse momentum observed experimentally in 
order to describe nuclear collisions. The question arises 
whether one can derive the parameters of the statisti- 
cal model of nucleus-nucleus collisions from the parame- 
ters of pp collisions with the geometry and the dynamics 
of parton rescattering as input. In this paper we have 
demonstrated that this is possible. 

Let us review our assumptions. The derivation of the 
average path length of a parton in the target nucleus uses 
the number of participants as weight factor. This num- 
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ber depends weakly on the inelastic pp cross section and 
therefore induces small differences in the centrality de- 
pendences of the average profile functions of Au-Au and 
Pb-Pb collisions at RHIC and LHC energies; cf. Fig. 1. 
Only a minor variation is observed in experiment. Per- 
haps here an improvement is possible. 

We consider the intrinsic momenta of the gluons to be 
described by the LCP distribution in order to average 
the dipole cross section and the gluon structure function. 
The LCP distribution comes from a statistical model for 
the nonequilibrium distribution of the gluons. It is an 
assumption that the intrinsic gluons have the same trans- 
verse momenta as the liberated gluons. 

For the mutual boosting estimate the shadowing effect 
is taken into account very crudely. Furthermore, it is 
subtle that these two effects are counteracting each other 
leading to the same result as the crude calculation from 
parton multiple scattering. This is not expected at all 
energies and should be explored in more detail. 

The simple features of the presented calculation are 
the most convincing: The multiplicity at pseudorapidity 
77 = is roughly proportional to the product ^^j.) 
[T]. Therefore, the geometrical dependence of trans- 
verse momentum broadening is clearly visible in the be- 
havior of the rapidity distributions. The same depen- 
dence of the rapidity distributions on centrality in Au- 



Au at RHIC and Pb-Pb at the LHC confirms this find- 
ing. There is a clear increase of the magnitude of the 
rapidity distribution at = for the different ener- 
gies. At LHC {y/SNN = 2.76 TeV) one obtains roughly 
a twice larger value for all centralities than at RHIC 
{^/snn — 0.2 TeV). This factor arises from the twice 
as large intrinsic mean (p^) in pp collisions (cf. Table 1) 
and from the twice as large momentum broadening Ap^ 
from multiple scattering reflected in the energy depen- 
dence of the transport parameter {crp]_). 

So our presentation has been able to separate the ge- 
ometrical from the dynamical effects. Very recently we 
received a preprint |20j simulating rapidity distribution 
at 77 = including event-to-event fluctuations. On aver- 
age, the above authors reproduce the same geometrical 
picture and saturation, including a microscopic picture 
of the initial gluon plasma fields. 

ACKNOWLEDGMENTS 

The work of B.Z.K. was supported by the Al- 
liance Program of the Helmholtz Association, contract 
HA216/EMMI "Extremes of Density and Temperature: 
Cosmic Matter in the Laboratory", and by Fondecyt 
(Chile) grant 1090291. 



[1] H . Pirner and K. Reygers, Ph ys.Rev. D86, 034005 
(2012), arXiv:1106.5486 [hep-ph] 



J. Bjorken, Phys.Rev, 
P. F. Kolb and U. 



D27, 140 (1983) 
W. Heinz, (2003), 



arXivinucl- 



1^ 



th/0305084 [nucl-th] 
S. Domdey, D. Grunewald, B. Kopeliovich, and 
H. Pirne r, Nucl. Phys. A825, 200 (2009) , arXiv:0 812. 2838 
[hep-ph] I 

L. D. McLerran and R. Venugopalan, 'Phys.Rev. D49 



"ITT 



2233 (1994), ar Xiv:hep-p h/9309289 [hep-ph]; Ph ys.Rev. 
p49, 3352 (1994) , arXiv:hep-ph/9311205 [hep-ph]* 
M. Johnson, B. Kopeliovich, and A. Tarasov, Phys.Rev. 
C63, 035203 (2001), arXiv:hep-ph/0006326 [hep-ph] 
B. Z. KopeUovich, A. Schafer, and A. V. Tarasov, 
JPhys.Rev. D62, 054022^^(2000 J| |arXiv:hep-ph/9908245| 
Jhep-ph] 

[8] H. De Vries, C. De Jager, and C. De Vries, Atom. Data 

Nucl.Data Tabl. 36, 495 (1987). 
[9] M. L. Miller, K. Reygers, S. J. Sanders, and P. Stein- 
berg, Ann.Rev.Nucl.Part.Sci. 57, 205 (2007)||arXiv:nucL] 



r 



10 



_ex/0701025 [nucl-ex]J" 
_S. Adler et al. (PH ENIX 
C71, 034908 (20057, 



Collaboration) , 



""arXiv: nucl-ex /04090 1 5 

K. Aamodt et al. (ALICE Collaboration), Phys.Rev. Lett. 



[Phys.Rev 

]nucl-exj 



106, 032301 (2011)HarXiv:1012.1657 [nucl-ex]| 



3tt] 



[11] 



TI2T 



B. Kopeliovich, I. Potashnikova, and I. Schmidt, 
Phys.Rev. C81, 035204 (2010), arXiv: 1001. 4281 [hep-, 

m 



B. 

and 



Kopeliovich, 
I. Schmidt, 



I. 



Potashnikova, 



[13] 



'arXiv: 1007. 1913 [hep-ph] 
B. Kopeliovich, I. Potashnikova, 

CSS, 



H. Pirner, 

iPhys.Lett. B697, 333 (2011) 



H. Pirner, and 
014912 (2011)1 



I. Schmidt, Phys.Rev 
[arXiv: 1008.4272 [hep -ph] 
B. Kopeliovich, I. Potashnikova, and I. Schmidt, 
'Nucl.Phys. A864, 203 (2011), arXiv:1012.5648 [hep-ph]J 
A. Martin, W. Stirhng, R. Thorne, and G. Watt, 
Eur.Phys.J. C63, 189 (2009), arXiv:0901.0002 [hep-ph]t 
J4. Cluck, E. Reya, and A. Vogt, Eur.P hys.J. C5, 46l] 
(1998), arXiv:hep-ph/9806404 [hep-ph], 
A. Dumitru, D. E. Kharzeev, E. M. Levin, and Y. Nara, 
Phys.Rev . C&b, 044920 (2012)1 | arXiv:1111.3031 [hep-j 

A. Shosh i, F. Steffen, H. G. Dosch, and H. Pirner, 
[Phys.Rev. D66, 094019 (2002)| |arXiv:hep-ph/0207287| 
[hep-ph( \ 



S. Albino, B. Kniehl, a nd C. Kram er, [NucLPhys. B725^ 
181 (2005), |arXiv:hep-ph/ 0502188 [hei3-ph[' 
B. Schenke, P. Tribedy, and R. Venugopalan, (2012), 
[arXiv: 1206. 6805 [hep-ph]] 



